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Abstract: 

The reduced velocity correlation function for fragments from the reaction Fe + Au at 100 
A MeV bombarding energy is investigated using the dynamical-statistical approach QMD+SMM 
and compared to experimental data to extract the Freeze-Out volume assuming simultaneous 
multifragmentation. It is shown that the data are consistent with a Freeze-Out volume corre- 
sponding to 0.1 —0.3 times normal nuclear matter density. The calculations show an additional 
correlation due to recoil effects from a heavy third fragment in an asymmetric break-up. This 
effect is present, but less pronounced in the data. 
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A topic of great interest in nuclear physics is the multifragmentation of heavy nuclei at mod- 
erate excitation energies. Such systems are typically formed in p-nucleus reactions in the GeV 
energy regime or in intermediate energy heavy ion collisions. In the experiments of the Purdue 
group first attempts were made to relate the inclusive results from p-nucleus fragmentation 
to the critical exponents of a phase transition, analogous to the liquid gas transition in con- 
densed matter physics. Within the last few years it has become possible to study such reactions 
in highly exclusive experiments in order to extract the critical properties. Critical properties 
can then be extracted investigating by, for example the moments of the mass distribution as 
proposed by Campi 0. 

Such a procedure, however can not provide the physical quantities that drive this phase 
transition. If one wants to extract quantities like a critical temperature and density one has 
to rely on dynamical models. Unfortunately there is presently no complete model available 
that describes the process of thermally driven multifragmentation in heavy ion collisions in 
a consistent approach. Current modeling involves a two-step dynamical-statistical process 
beginning with a pre equilibrium stage described in the earliest approaches by an Intranuclear 
Cascade Model and, more recently, by molecular dynamics pictures such as the Quantum 
Molecular Dynamics model, QMD |^ or by single particle models of the VUU/BUU type 
0. Following this stage there remains a distribution of nucleons and complex fragments which 
can themselves be highly excited and will undergo further statistical decay. This statistical 
decay process has been modeled in various codes as a sequential evaporation R, & M or as an 



explosive simultaneous multifragmentation |TT|, |12|, None of these models, however, 

includes the dynamics of the fragmentation. 

Recent investigations with the QMD model have shown that for central collisions of heavy 
nuclei at high energies, a rapid compression-decompression mechanism emerges in the initial 
stage of the reaction and leads to a direct multifragmentation process, where a highly excited 
heavy residue is no longer formed with high probability P, |l5l. These investigations suggest 
that the region to search for the thermally driven multifragmentation in heavy ion reactions is 
in central collisions at moderate bombarding energies (E/A ^ 100 A MeV) or in peripheral (or 
highly asymmetric) reactions at higher energies (E/A «i 1 A GeV) . Under these conditions the 
direct reaction leads to a highly excited source that then breaks up into many IMF's. 

Recently, attempts have been made by Sangster et al., and Lacey at al. to extract 



the emission time scale of the multifragmentation process by assuming a sequential decay 
mechanism with a fixed Freeze-Out volume. In contrast to earlier work performed by Kim et 
al. using the Koonin-Pratt formalism []19| they used a the classical three-body trajectory 



code MENEKA j^/ to simulate the sequential emission of IMF's from the surface of a spherical 
source characterized by a unique radius parameter. In this Letter we will show, that the same 
data can be explained by assuming a simultaneous multifragmentation from a Freeze-Out 
volume considerably larger than the ground state dimension. 

This will be done by using the QMD-I-SMM approach. In this case the Hamilton equations 
of all nucleons are integrated, which implies that correlations between all the existing nucleons 
and fragments are treated in all orders. This is especially important in the case when one 



1 



deals not only with two IMF's, but with many. In the following we present comparisons of the 
QMD+SMM fragment-fragment correlation function for the reaction Fe (100 A MeV) + Au 



for central collisions (b = - 6 fm) in to the data of Ref. After the break-up of the hot 



target remnant within the SMM model (after 300 fm/c reaction time) we neglect the nuclear 
forces and follow the Coulomb trajectories of the charged particles through an experimental 
aceptence filter. 

The QMD+SMM approach models the dynamical reaction and subsequent statistical decay 
of the excited pre fragments (details about the two models can be found in refs. |12|). 
Previous investigations with the QMD model ^ have shown that there are two different mech- 
anisms leading to the multifragmentation. One is related to the mechanical rupture of the 
system whenever compressional effects are important; the other produces fragments thermally 
from an equilibrated source. This thermal multifragmentation has so far not been described in 
a microscopic model like QMD ||^. In the first comparison to inclusive IMF data it was shown 



that a two-step model was necessary to reproduce the experimental angular distributions [21 



This two-step model involved the calculation of initial kinematics and excitation energy of all 
pre fragments with the QMD model followed by a subsequent deexcitation calculation utilizing 



the Statistical Multifragmentation Model (SMM) of Botvina et al. [jT2|. The input for the SMM 
stage of the reaction, i.e., the mass and the excitation energy of the fragments, are consistently 
determined within the QMD approach. 

The SMM model describes the multifragmentation of highly excited nuclei based on the 
statistical approach and a liquid-drop description of hot fragments. It is assumed that the 
excited primordial fragments break up into an assembly of nucleons and fragments. All these 
decay products are described as Boltzmann particles in a Freeze-Out volume = Vo(l + k), 
where k is a model parameter and Vq is the volume of the system corresponding to normal 
nuclear matter density. Since all the produced fragments are excited (only particles with A < 
4 are considered as elementary particles) the final deexcitation is treated as a Fermi break up 
for lighter fragments ( A < 16) and via an evaporation of nucleons and clusters up to ^^O for 
heavier fragments (for details see Ref. |]12|). In the following we will vary the volume parameter 
K from 1 to 15 in order to extract the Freeze-Out volume by comparing the reduced velocity 
correlation functions for each value of k with the experimental data. 

Figure |T] shows a comparison between the experimental data |2T] and the QMD+SMM 
calculations using k = 2, 5 and 10 for two different projections of the double differential cross 
section, d?a / dZ / dVL. Both the data and the calculations are for central collisions only (for details 
see Ref. [0). The upper part of figure |1| shows the charge yield distribution at a laboratory 
angle of 72° (±13°), while the lower part shows the angular distribution of fragments with 
Z=10. In both cases it can be seen that the variation of the volume parameter k does not 
influence these semi-exclusive observables and all calculations agree reasonably well with the 
data. However, the calculated angular distributions are still somewhat too steep and under 
predict the data at backward angles. 

The reduced velocity correlation functions shown in figure |^ have been obtained by taking the 
ratio of the correlated reduced velocity distribution and a background distribution obtained by 
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event mixing between correlated IMF's in a single event divided by the same quantities obtained 
from two IMF's from different physical events, thereby eliminating final state correlations. In all 
cases the true (Ytrue{vred)) and background distributions (Yhack{vred)) were developed separately 
for two classes of events: 1) two IMF's detected on opposite sides of the beam and, 2) two IMF's 
detected on the same side of the beam. In the areas where these two correlation functions 
overlap in relative velocity it was found that the correlations obtained were effectively identical 
and so in this paper the results have, in all cases, been combined into single distributions. The 
two fragment correlation function is then calculated according to 

Fig. shows the experimental correlation function (symbols) for the three different charge 
product bins Zi • Z2 = 25 — 64 , 65 — 129 and 130 — 250 (the smallest detected charge is Z=5 
in all cases) compared to the calculated results (lines) for different volume parameters k in the 
SMM model (recall that V = Vo(l + k)). The curves show fits to the theoretical correlation 
functions using the fitting function 

l+R{Vred)=a ^ j^ . (2) 

1 + e c 

It can be seen that, in the limit of simultaneous fragment emission, the size of the Coulomb 
hole can be explained within the QMD+SMM approach if a volume parameter between k = 2 
and 10 is used; both the results with n = 1 and k = 15 clearly disagree with the data 
for the heavier fragments. This means that these correlation functions are consistent with a 
simultaneous break-up with a Freeze-Out density of ^ = 0.1 — 0.3 go- For the lightest fragments 
the smallest volume parameter seems to describe the data best, but here the sensitivity to the 
Freeze-Out volume is less pronounced. This may indicate that the fragments do not come from 
a single Freeze-Out volume, but that the smaller fragments are emitted earlier from a smaller 
volume. 

These results are in good agreement with previous investigations of the same data using the 
sequential three-body trajectory code MENEKA [|l^. In this case an emission time of order 
500 fm/c was found for a fixed freeze-out volume equivalent to k = 1. These time scales are 
comparable to the typical time a fragment needs to traverse a Freeze-Out volume of roughly 5 
times the nuclear volume. 

Another feature that can be observed in the calculated correlation functions in figure ^ is 
the pronounced enhancement at a Vred ~ 15 — 20. This additional correlation is absent in the 
present light fragment data and is not seen in the data of the MSU group except for very 



peripheral collisions p2[. Our calculations also show that this effect is more pronounced for 
larger impact parameters and increases (see figure |^) with decreasing SMM volume parameter 
(a similar behavior has recently been found in the Berlin model when the excitation energy 



is decreased). A similar correlation was recently observed for the Au + Au reaction at 150 
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A MeV bombarding energy [Q. In this case the enhancement has its origin in the collective 
flow of the fragments. In our case, however we are dealing with a more asymmetric reaction 
at lower energies where the collective flow vanishes due to the compensation of the repulsive 
and attractive nuclear forces. However, the Coulomb repulsion due to a heavy third fragment 
can produce an additional correlation in the fragment-fragment correlation function. In table |l| 
we show the average charge of the largest fragment as well as the average charge asymmetry 



< AZ >= ^((Zi - Z2Y + (Zi - Z'iY + (Z2 - Z3)2)/3 and the average multiphcity of IMF's 
with Z=4-20 for values of k from 1 to 15 over the impact parameter interval 0-6 fm in the reaction 
Fe (100 A MeV) + Au. It can be seen that the charge distribution gets very asymmetric with 
decreasing k and the average IMF multiplicity decreases by a factor of two (the same holds 
also when the impact parameter is increased for a fixed volume parameter). This is a strong 
indication that the enhancement in the calculated correlation functions in figure |] results from 
a large charge asymmetry. 

The fact that the experimental data in ref. [|I^ does not show such an enhanced peak 
indicates that the break up pattern in nature is more symmetric than described within the 
SMM model. In order to clarify this point we show in figure ^ the experimental data for events 
with three detected IMF's. For this subset of the data, only the two lightest fragments are 
used to generate the correlation functions; the heaviest fragment is simply tagged as Zmax- The 
curves in Fig. ^ have been generated based on the following Zmax selection criteria: a) the 
numerator in eq. |l] contains only pairs from events with Zmax > 18 (note however that this 
refers to the largest detected fragment) (dotted line), while no requirement has been applied 
to the fragments in the background pairs (the denominator in eq. 0); b) both the correlated 
and the background pairs are from events with detected Zmax > 18 (dashed line); c) no cuts 
on Zmax (full line). 

One clearly observes that, in accordance to the QMD+SMM results, a large peak is found for 
the trigger condition a) due to the recoil of a heavy third fragment which is not accounted for in 
the background. This peaks goes away when the background fragments are taken from events 
with a similar charge asymmetry. The full line shows again that the enhancement virtually 
disappears when no cuts are made on Zmax- 

We have shown that the fragment-fragment reduced velocity correlation functions can be 
explained within the two-step model QMD+SMM, where the initial, dynamical step of the re- 
action, modeled using the microscopic QMD approach, is followed by a simultaneous statistical 
multifragmentation. The correlation functions for the system Fe (100 A MeV, b=0-6fm) + Au 
are consistent with a simultaneous break-up with a Freeze-Out density off) = 0.1 — 0.3^o- 
These findings are consistent with previous investigations assuming a sequential emission of the 
fragments from a fixed volume source. The typical time seal for this process was found to be less 
than 500 fm/c. We note that these time scales are comparable to the typical time a fragment 
needs to traverse a Freeze-Out volume of roughly 5 times the nuclear volume. The calculations 
show a additional correlation due to the coulomb repulsion from a heavy third fragment in an 
asymmetric decay. This effect is less pronounced in the experimental data, unless one explicitly 
triggers on pairs from events with a large charge asymmetry. 
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Figure Captions: 



Figure 1: 

Projections of the semi-exclusive triple differential cross section (Pa/dZdVL for the reaction Fe 



1^100 A MeV) + Au. The symbols show the data for the fragment charge distributions 



at "dhab = 72° (upper panel) and the fragment angular distributions for Z = 10 (lower panel). 
The histograms show the calculations with the QMD+SMM (b = - 6 fm) model for different 
volume parameters, indicated. 



Figure 2: 

Mixed fragment reduced velocity correlation functions for three different constraints on the 
Coulomb product Zi ■ Z^- The symbols show the data of ref. [|TB|, while the curves are fits to 
equation |^ obtained with the QMD+SMM model. The different lines show calculations done 
with the volume parameter k = \ (solid line), k = 1 (variable length dashed line), k = 5 
(dashed line), /t = 10 (dash-dotted line) and k = 15 (dotted line). 



Figure 3: 

Experimental correlation function for the lightest two fragments in events with three detected 
IMF's. The full line shows the fit to the correlation function obtained from all events while for 
the dotted line includes only events in which the largest of the three detected fragments has 
a charge greater than 17. The dashed line shows the fit when both the true and background 
pairs are from events with Z^aax > 18 (dashed line). 
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Tables: 



K 




< AZ > 


< MULz=A-2Q > 


1 


26.5 


12.3 


2.2 


2 


22.7 


10.5 


2.7 


5 


18.8 


8.2 


3.4 


10 


17.9 


7.5 


3.6 


15 


16.0 


6.2 


4.0 



Table 1: 

The average charge of the largest fragment as well as the average charge asymmetry < AZ >= 
^((Zi - ^2)2 + (Zi - Z^Y + (^2 - ^3)2)/3 and the average multiplicity of IMF's with Z=4-20 
for the reaction Fe (100 A MeV, b=0-6 fm) + Au and values of k from 1 to 15. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9312008vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9312008vl 



This figure "fig3-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9312008vl 
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